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A thorough characterization of many samples of the filled skutterudite compound PrOs4Sbi2 is 
provided. We find that the double superconducting transition in the specific heat (T c i ~ f .89 K 
and T C 2 ~ 1.72 K) tends to appear in samples with a large residual resistivity ratio, large specific 
heat jump at the superconducting transition and with the highest absolute value of the specific 
heat above T c \. However, we present evidence which casts doubt on the intrinsic nature of the 
double superconducting transition. The ratio of the two specific heat jumps acIt^) snows a wide 
range of values on crystals from different batches but also within the same batch. This ratio was 
strongly reduced by polishing a sample down to 120^im. Remarkably, three samples exhibit a 
single sharp transition of ~15 mK in width at T c ~ 1.7 K. The normalized specific heat jump 
C — C normal I C normal) T of two of them is higher than ~ 32% so larger than the sum of the two 
specific heat jumps when a double transition exists. As an evidence of better quality, the slope 
in the transition is at least two time steeper. We discuss the origins of the double transition; in 
particular we consider, based on X-ray diffraction results, a scenario involving Pr- vacancies. The 
superconducting phase diagram under magnetic field of a sample with a single transition is fitted 
with a two-band model taking into account the good values for the gap as deduced from thermal 
conductivity measurements. 

PACS numbers: 65.40.Ba,71.27.+a,74.25.Dw,74.70.Tx,74.20.Rp 
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I. INTRODUCTION 



Since the discovery of the first Pr-based heavy fermion 
superconductor PrOs 4 Sbi2 (T c ~ 1.85K) by Bauer et 
ali, this system has attracted much attention with par- 
ticular emphasis on the possible unconventional nature 
of superconductivity. A significant piece of the evidence 
for unconventional superconductivity is the double super- 
conducting transition seen in specific heat first reported 
in 2003 by Vollmer et al3 and Maple et al4. Ever since, 
a plethora of publications have dealt with its observa- 
tion and with possible theories. This double transition 
has since been observed by specific heat measurements 
by many groups (from Japan 4 , USA 3 ' 5 , Germany^ 7 -^, 
France 1 ^.) and by thermal expansion^ 1 - with samples from 
different origins and even in La doped or Ru substituted 
samples 5 -*^. So whatever the origin, the double transi- 
tion is a robust property of this compound and we will 



see that it appears in good samples. 

However susceptibility measurements on a sample with 
a very clear double superconducting transition 10 induced 
the first doubts about its microscopic origin. Indeed even 
in this good sample, the diamagnetism is not perfect 
at T c i, the highest superconducting transition temper- 
ature, and a second step at T c2 , the lowest one, is vis- 
ible. This caused us to look more closely at other ev- 
idence for a microscopic origin, including the magnetic 
field, and pressure dependences of the double transition. 
The phase diagram under magnetic field didn't bring ev- 
idence for its intrinsic nature since the field dependence 
of T C 2 is completely similar to that of T c \. The shape 
of H C 2{T) can be quantitatively understood with a two- 
band models and indeed the same model simply scaled 
with T c has been used to fit both lines. The behavior of 
T c i ~ T C 2 under pressure is also not conclusive^ 3 -. At low 
pressure, the slope dT c \/dP is at least 20% smaller than 
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dT C 2/dP. However, above 1 GPa the behavior of the two 
transitions is similar, with T cl — T c2 stabilizing around 
200 mK. These results do not rule out an intrinsic ori- 
gin but certainly provide no supporting evidence towards 
it, (contrary to the well documented case of UPt3 where 
the different field and pressure dependences of the two 
transitions were decisive results). 

We report here on a study of the nature of the dou- 
ble superconducting transition of PrOs4Sbi2 . Our main 
purpose is to clarify whether the double superconducting 
transition which appears in specific heat is intrinsic, like 
it is now admitted for UPt 3 or extrinsic, due to sample 
inhomogeneities, as shown for URu2Si2 14 and high-T c su- 
perconductor YBa 2 Cu3 0^ 15 , We discuss this point from 
systematic characterizations by resistivity, specific heat 
and susceptibility measurements. From general charac- 
terizations, we conclude that the double transition ap- 
pears in good samples. But an extensive study partic- 
ularly of many small samples provides evidence which 
brings strong doubts about its microscopic origin; the 
most convincing one is the existence of three samples 
with a single sharp superconducting transition with a T c 
matching T C 2- Then we present the single crystal x-ray 
diffraction results. Finally we show the first measurement 
of the phase diagram under magnetic field for a sample 
with a single sharp superconducting transition, and a fit 
of the upper critical field with a two band model. 

We would like first to introduce some criteria we will 
use in this paper. Figure [1] provides the main criteria de- 
pending on the kind of superconducting transition (dou- 
ble, single and broad, or single and sharp), i.e. the su- 
perconducting transition T c , T c \ and T c2 obtained on the 
onset, the specific heat jumps A(C/T)i at T c ±, A(C/T) 2 
at T C 2 and A(C/T) for the whole jump, the width of the 
transition A(T cl ) at T cl , A(T c2 ) at T c2 or A(T C ) when 
a single transition appears, and the slope in the transi- 
tion a. This last criterion appears more relevant than 
the simple width of the transition when we compare one 
of the jumps in the double transition with a single sharp 
jump, which is probably two times higher. The residual 
resistivity ratio RRR is always measured between 300 K 
and 2 K. 



II. GENERAL CHARACTERIZATIONS 

The crystals of PrOs4Sbi2 were grown by the Sb flux 
metho d 16 i 17 ' 18 by 3 separate groups (P. Canfield, H. Sug- 
awara and G. Lapertot, respectively labelled C, S and L). 
The first number indicates the batch. When there were 2 
different samples from the same batch we add an extra in- 
dex. The crystals from P. Canfield and G. Lapertot were 
separated from flux by a hot spinning process. Remaining 
flux droplets were dissolved in hydrochloric/ water solu- 
tion. The samples used in ref i 10 ) 19 ; 20 ^ are respectively 
SI, L1-1A, S3, Ll-6. 

The specific heat (C) measurements were performed 
in a 3 He calorimeter either by a quasi-adiabatic method 
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FIG. 1: (color online) Definition of the criteria we use in this 
paper: (a) for the samples with a double superconducting 
transition, (b) with a sharp single transition and (c) with a 
broad single transition. 
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FIG. 2: (color online) Specific heat as C/T versus T for several 
samples of PrOs4Sbi2 from different batches. The samples' 
labels follow the decreasing order of the specific heat at 2 K. 
The absolute value of C clearly depends on the sample, with 
a value on the Schottky anomaly at 2 K included between 
2.06 J/K 2 .mol and 2.94 J/K 2 .mol. The insert is a zoom on 
the superconducting jump for the samples with the highest 
specific heat. A double superconducting transition appears in 
some samples while the other samples exhibit a broad single 
transition. 
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with a Au/Fe-Au thermocouple controlled by a super- 
conducting quantum interference device (SQUID) or by 
a heat-pulse relaxation technique using a Physical Prop- 
erty Measurement System (PPMS) from Quantum De- 
sign. 

Figure [2] presents the specific heat of several samples of 
PrOs4Sbi2. Their absolute values at 2 K differ strongly, 
varying from 2.06 J/K 2 .mol to 2.94 J/K 2 .mol and can 
vary on the samples in the same batch (see samples L2-1 
and L2-2). Other published results report absolute val- 
ues at 2 K between 1.3 J/K 2 .mol 3 and 3.2 J/K 2 .mol 8 . 
This discrepancy cannot be explained only by the pres- 
ence of trapped Sb flux. The specific heat broad peak at 
2 K is mainly due to the Schottky term due to the pres- 
ence at ~ 8K of the first crystalline electric field (CEF) 
excited level of the Ai 2 Pr states, I4 , above the singlet 
ground state Ti . As we point out to possible Pr- vacancies 
in the structure (see section III.B), we propose that the 
samples dependence of the specific heat is partly due to 
these Pr-vacancies. Quantizing the necessary percentage 
of vacancies to fit the curves is a difficult task. Indeed, 
some distortions may appear with Pr-vacancies and they 
may change locally the CEF, which broaden the Schot- 
tky term. Moreover the interactions between Pr ions are 
quite strong as shown by the dispersion of the inelatic 
spectrum in wavevectors^; so the consequence of the Pr- 
vacancies cannot be reproduced by a pure local model. 
Roughly a maximum of 10 % of Pr-vacancies is required 
to fit the specific heat of sample L3-1 to Ll-1. No upturn 
due to the nuclear Schottky anomaly except a small fea- 
ture for the samples of batch CI, was observed, at least 
above 0.4 K. We also note that no anomaly was detected 
at 0.6 K, temperature at which several experiments re- 
port a change of behavior—. 

The double superconducting transitions for several 
samples are shown in the zoom of Fig. [5] (see insert) . The 
shape of the double transition of samples from different 
origins (Ll-1 and Sl-2) is similar but is quite different 
from sample L3-1 with similar T c s and even larger devia- 
tions are found from sample L2-2 both in absolute value 
of C and in value of T C 2 . The sample dependence of the 
shape of the double transition will be discussed in detail 
in section UTTl 

We discuss now the link between the appearance of 
the double transition and sample quality. Figure [3] shows 
the superconducting transition temperatures determined 
by specific heat and resistivity measurements versus the 
residual resistivity ratio. In both cases T c was obtained 
from the onset criterion. When the specific heat and the 
resistivity were not measured on the same sample (which 
was often the case), the RRR values taken are an aver- 
age for the batch. As the spread of the RRR is large 
(for instance the RRR of batch LI is included between 
17 and 33.8), only a general tendency can be extracted. 
However it is quite clear that T c \ increases slightly with 
the RRR (its minimum and maximum values are respec- 
tively 1.805 K and 1.897 K) whereas T C 2 is strongly sam- 
ple dependent. All the smaller T C 2 at about 1.5 K were 
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FIG. 3: Superconducting transition temperature versus the 
RRR between 300 K and 2 K determined by specific heat 
(full square and circle) or resistivity (cross) measurements. 
The open circles and squares show the T c s determined by C 
and p measurements of the samples with a single sharp su- 
perconducting transition Ll-5 and L1-1A. Definition of the 
T c s is shown in the Fig [T] The large circle includes all the 
samples with a clear double superconducting transition (ex- 
cluding batch L2) and sample L1-1A. The general tendency 
is that the double transition appears in the samples with the 
highest RRR. 



reported from samples belonging to batch L2. The large 
circle indicates all the samples with a clear double transi- 
tion (excluding sample L2-3) and also contains a sample 
with a single sharp transition L1-1A. The double tran- 
sition seems to appear in the samples with the largest 
RRR. We notice that the samples with a single sharp 
transition have a T c similar to T C 2 and a relatively high 
but not the highest RRR. They will be discussed further 
in section III. 



Figure 4(a) shows the absolute value of C on the Schot- 
tky anomaly at 2 K against the RRR. The highest the 
RRR is, the largest the specific heat at 2 K is. Fig- 
ure |4(b)| presents the total jump of the specific heat at the 
superconducting transition (A(C/T)), which is probably 
the best criterion of the quality of the samples, against 
the value of C at 2 K. Again, the highest the specific heat 
jump is, the largest the specific heat at 2 K is. Moreover 
as indicated by the large circles on Fig. 4(a) and |4(b)| 
which embodies all the samples with a double transition, 
the double transition seems to be a feature of the samples 
which meet all these criteria. Samples with a broad sin- 
gle superconducting transition, which always covers the 
temperature range of the double transition, are clearly of 
worse quality than the double transition samples. 

However, we will present in the next section some 
characterisations which provide several pieces of evidence 
against an intrinsic origin of the double transition, spe- 
cially those carried out on very small samples (with a 
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FIG. 4: (a) Absolute value of the specific heat C/T on the 
Schottky anomaly at 2 K versus the RRR. The highest the 
RRR is, the larger the specific heat at 2 K is. (b) Entire 
specific heat jump A(C/T) (defined figure [TJ versus the spe- 
cific heat C/T at 2 K. The two large circles embody the sam- 
ples with a clear double superconducting transition (excluding 
batch L2) : the samples with a double transition are charac- 
terized by the highest C at 2 K and the highest specific heat 
jump at the superconducting transition. 



typical size of 100 fim and lighter than 0.1 nig). 



III. DOUBLE SUPERCONDUCTING 
TRANSITION 

A. Characterization results 

The first doubts about the intrinsic nature of the dou- 
ble superconducting transition came from the results of 
susceptibility measurements. Indeed no published result 
shows perfect diamagnetism at T c \ : the ac susceptibility 
always exhibits a very broad superconducting transition 
or a double step matching with T cl and T c2 of the specific 
heat result s 7 ' 10 . All the samples with a double transition 
we have tested by ac susceptibility (x), i.e. batches LI, 
L2, SI, also exhibit a double step matching the specific 
heat jumps. We report here other stronger doubts, par- 
ticularly based on the existence of samples with a single 
sharp superconducting transition. We will compare their 
quality with those of samples with a double transition. 

Table [J provides a general view and the main features 
of all the samples with a double transition or a single 
sharp superconducting transition we have tested. All of 
them except S1-1A and L1-1A are as grown (unpolished) 
samples. All the samples not reported in table U exhibit 
a single transition with a width larger than the temper- 
ature range of the double transition as shown in Fig. [5] 
The criteria for T ci , A(T ci ), A(C/T)i, a Tci are defined in 
the Fig. [0 

We first notice that the double transition appears in 
5 different batches from 3 origins with similar shape for 
the biggest samples of the batches LI, SI and CI. But 
further investigations, specially on small samples (with 
a typical size of 100/im), show a wide range of the ratio 
A(C/T) 1 /A(C/T) 2 from 0.2 to 3.36. It even differs in 
the same batch, from 0.2 to 0.7 or from 0.25 to 3.36 in, 
respectively, batches LI and L2. If we suppose that one 
of the superconducting transitions is more sensitive to 
the quality of the sample than the other, other proper- 
ties should also be affected. However it clearly appears 
that we cannot connect the value of the ratio with other 
criteria, as the value of T c2 or as the width of the tran- 
sition at T c2 , At c2 , by comparing, respectively samples 
L2-2 and L2-3 or samples L3-1, Ll-1 and Ll-6. 

Before discussing further the characterizations, we fo- 
cus on our most remarkable finding which is shown in 
Fig. O we measured the specific heat of three samples 
Ll-5, Ll-7 and Ll-lA with a single sharp superconduct- 
ing transition (Ll-lA has still a tiny jump at 1.80K as 
well as a small step in resistivity at 1.85K 19 ). The specific 
heat of a sample with a "usual" double transition (Ll-1) 
is also presented. We notice the existence of sample L2-3 
with a clear double transition but with a quite reduced 
T c2 . Ll-5 and Ll-7 are very small as-grown platelets 
with well-developed faces and with a thickness of about 
50/im. Ll-lA has been obtained by polishing a large 
cube (1 mm) of the sample Ll-1 so that the thickness 
was reduced down to 45 fim^. Their critical tempera- 
ture T c is 1.733 K, 1.680 K and 1.745 K, so matching T c2 
of the samples with a double transition. One of these 
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sample 


shape _ size/weight _ measurement 


Tcl(K) 


Tc2(K) 


A(C/T)l 


A(Tcl)(mK) 


A(T c2 )(mK) 


A(C/T) 






A(C/T) 2 


C/T normal 


Double Superconducting Transition 


Sl-2 
Ll-1 
Ll-3 
Ll-4 
Ll-6 
L3-1 
L3-2 


AofC. „ 1.97mg 
AofC. - ~10mg 
c. „ 200^m 
b. „ 500/xm 

p. „ 50*150*150pm 3 „ AC 
AofC. „ 1.25mg 
b. „ 200/im 


1.884 
1.890 
1.868 
1.889 
1.85 
1.891 
1.897 


1.721 
1.737 
1.74 
1.732 
1.73 
1.761 
1.76 


1.23 
0.70 
0.29 
0.20 
0.54 
2.8 
1 


61 

61 
43 
64 
36 
46 
60 


23 
35 
40 
35 
39 
21 
40 


30.6% 
32.3% 
30.7% 

28% 
>27.2% 
27.3% 

32% 


4.1 
4.23 
4.7 
5.2 
3.7 
3.7 
4.1 


1.98 
1.58 
1.3 
0.79 
2.0 
3.4 
2.14 


Sharp Single Superconducting Transition 


Ll-5 

Ll-7 
L1-1A 


p. „ 160*200*40/im 3 „ PPMS 
„ AC 

p. „ 50*150*150Aim 3 „ AC 
p.„ 45*150*200/im 3 /~0.1mg 





1.800 


1.733 

1.680 
1.745 





~0 (<0.09) 





55 


17 

21 
27 


33%< <44% 
>28% 
>22% 
36% 


16±4 
16.5±2 
8.1 
10±2 








Double Superconducting Transition with a weaker T C 2 


L2-2 
L2-3 
Cl-2 


b. „ 9.7mg 

b. „ 150/im 

c. „ 200/im 


1.844 
1.85 
1.877 


1.53 
1.535 
1.685 


3.36 
0.25 
1.17 


150 
75 
95 


105 
47 
35 


16% 
28% 
27% 


0.3 
4.6 
2.2 


1.05 
1.75 
1.4 


Effect of polishing : change of the relative height of the two specific heat jumps 


Sl-1 
S1-1A 


p.„ 300/im „ AC 
p. „ 120pm „ AC 


~1.88 
~1.88 


~1.71 
~1.71 


1.9 
~1.0 


100 
180±20 


25 
20 


>28% 
>28.1% 


4.2 
7.74 


1.5 
1.0±0.3 



TABLE I: Main features of the specific heat of the samples with a double superconducting transition or a single sharp tran- 
sition. To get S1-1A and L1-1A, Sl-1 and one piece of Ll-1 were respectively polished. Abbreviations are: c.=cube; b.=bar, 
AofC.=aggregate of cubes, p.= platelet; AC=ac calorimetry measurements. The slope at T C i, ar ci is taken on the specific heat 
jump normalized to its normal value. The criteria are indicated in the Fig. [T] 

A large range of ratio of the two specific heat jumps appears. Three samples with a single sharp transition exist: Ll-5, Ll-7 and 
L1-1A. The superconducting transition of sample Ll-5 is the sharpest. For samples Ll-5 and L1-1A and relatively to the specific 
heat in the normal state, the specific heat jump is higher than the sum of the two specific heat jump in the double transition 
samples. The slope in the superconducting transition at T C 2, ar a2 , increases with decreasing ratio A(C/T)i/A(C/T)2, reaching 
the highest values for the 3 samples with a single transition (the transition is even three times steeper for the single transition 
sample Ll-5). 



samples Ll-5 was further characterized by single crystal 
x-ray diffraction, ac susceptibility and resistivity mea- 
surements. We confirmed the composition of the sample 
(see section III.B). For the susceptibility measurement a 
tiny susceptometer in glass fiber was built to get a good 
filling ratio. The secondary coil has an external diame- 
ter of 500 pm and hole of 300 /im. 440 turns of 14 pm 
diameter copper wires were wound on each part on the 
secondary coil. The frequency of the exciting magnetic 
field of 0.36 mT is about 375 Hz. Ac susceptibility and 
resistivity results are visible in Fig. [6j From these results, 
we checked that this sample has a single superconducting 
transition. Indeed, the susceptibility exhibits no sign of 
superconductivity above T c . All T c (p, Xt C) are consis- 
tent. 

Now we compare the quality of these single transition 
samples with the samples with a double transition. As 
the quantitative value of the specific heat of most of the 
samples is indeterminable (because many samples are too 
small or were measured by ac calorimetry), we have sub- 
tracted the normal part obtained under magnetic field 
up to 1 T then normalized to this value. We are only 
able to consider semi-quantitative values because, as we 



discussed in section II, the absolute value of the specific 
heat depends strongly on the sample. The ac calorime- 
try results of samples Ll-7 and Ll-6 provide only a bot- 
tom limit of the normalized specific heat jump as the 
background was not subtracted. As for the sample Ll-5, 
a semi-quantitative measurement of C (its mass, about 
0.1 mg, was too small to be measured precisely) was per- 
formed with the PPMS (Cf. figure [5]). The most impor- 
tant feature is certainly that the normalized specific heat 
jump, AC/C)t c , is at least as high as the entire tran- 
sition of the double transition samples. Indeed even if 
the smallness of the sample Ll-5 implies a noisy specific 
heat measurements by the PPMS (shown in Fig [5]), the 
specific heat jump at T c was found to be 33-44% of the 
normal state specific heat, a value larger than the entire 
jump, AC/C)t c , in any sample showing a double transi- 
tion (maximum of AC/C)t c =32.2%). The conclusion for 
Ll-1 A is similar with a AC/C)t c of 36%. The specific 
heat of Ll-7 was only measured by ac method without 
subtracting the background; the specific heat jump is un- 
derestimated. 

As for the width of the transition, sample Ll-5 has 
the sharpest ever measured (16 mK in p, 17 mK in C 
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FIG. 5: (color online) Zoom on the superconducting tran- 
sition of the specific heat normalized to its normal value of 
selected samples of PrOs4Sbi2. Samples Ll-1 and L2-3 ex- 
hibit a double transition, with the most common shape for 
Ll-1 and with a quite different T C 2 for L2-3. Our remarkable 
finding is the existence of three samples Ll-5, Ll-7 and Ll- 
1A with a single sharp transition (with a tiny remaining jump 
at 1.8 K for L1-1A). It clearly questions the intrinsic nature 
of the double transition. Two single transition samples Ll- 
5 and Ll-7 were measured by an ac method. Their specific 
heat jumps are underestimated. The specific heat of sample 
Ll-5 was also determined semi-quantitatively by relaxation 
method (Cf. figure (SJ. All features are reported in table U 
and comparison of the quality of the samples is provided in 
the text. 



and 35 mK in x)- One remarkable feature is that the 
transition is three times steeper (ar c three times larger) 
than the steepest transition in the samples with a double 
transition. The values obtained for a may be influenced 
to some extent by the range of values of C in the normal 
state, but insufficiently to change this semi-quantitative 
conclusion. Actually, this tendency is confirmed for all 
the samples of table [IJ the ratio A(C/T)i/A(C/T) 2 de- 
creases with increasing steepness of the transition at T C 2 
reaching zero for the three steepest transitions. This cri- 
terion as well as the specific heat jump AC/C)t c points 
to a higher quality of the samples with a single transition. 

We think that the transition at T C 2 is the intrinsic one. 
Indeed the transition at T c \ is always broader and less 
steep than the transition at T C 2 (Cf. table HJ. All the 
single sharp transition samples have T c lower than 1.75 K 
whereas even in the worst samples (see Fig. [3]) T cl is not 
smaller than 1.805 K. From this observation, we exclude 
that the single transition is the transition at T cl shifted 
to smaller temperature by some impurities effects. It 
would also be quite surprising that in these single transi- 
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FIG. 6: Ac susceptibility and resistivity of the sample Ll-5. 
The width of the transition in resistivity is 16 mK and 35 mK 
in susceptibility. T c is consistent for all measurements (C, p 
and x)- The RRR is 26. There is no trace of superconduc- 
tivity above T c . This sample has a single superconducting 
transition which happens to be the sharpest in all measure- 
ments. 



tion samples the transition at T c \ simply disappears due 
to worse quality of the samples; because the single tran- 
sition should also be broadened and because the thermal 
conductivity (k) measurements clearly point to a higher 
quality for sample with a single transition called L1-1A 
(k/T~ 70 ^W/K 2 .cm at 100 mK 19 ) than for a sample 
which exhibits a large superconducting transition with 
T c about 1.85 K (k/T~ 250 /iW/K 2 .cm at 100 mK 2 !). 

Finally only two points could leave some doubt. First, 
the RRR of samples Ll-5 and L1-1A are respectively 26 
and 30 whereas the largest value we have got (43) is in a 
double transition sample (Sl-1). Of course, the resistiv- 
ity is not a probe of the whole volume of the sample. For 
instance, sample Sl-1 exhibits a sharp superconducting 
transition in resistivity (p=0 at 1.82 K) but two steps in x 
matching with T c \ and T C 2- Moreover we cannot exclude 
that another parameter, such as remaining flux or Pr- 
vacancies, has the opposite effect on the RRR. Secondly, 
if we assume that the intrinsic T c is T C 2, we can compare 
T c of the single transition samples with T C 2 of the double 
transition samples. It appears that the transition in the 
single transition samples is shifted to smaller tempera- 
ture than the transition at T C 2 of the double transition 
samples of batch L3 (T c is smaller than T C 2 and the tran- 
sition in batch L3 is not much broader). One can argue 
that the parts of the sample with different T c may be af- 
fected independently by others impurities not related to 
the appearance of the double transition. This hypothesis 
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would also explain the existence of such samples as L2-2 
and L2-3 with a T c2 15% smaller than the highest one 
but with similar T c \. Finally, the rareness of the samples 
with a single sharp transition (to our knowledge, no other 
observation was reported) is certainly due to their small- 
ness and all the consecutive experimental difficulties. 




Sl-1 polished down to 120 mm 



1.4 1.5 1.6 1.7 1.8 1.9 2 2.1 2.2 

T(K) 

FIG. 7: AC specific heat versus temperature of samples Sl-1 
and S1-1A we got after polishing Sl-1 down to 120/im. The 
ratio of the two specific heat dramatically changes from 1.9 
to ~1.0. 



Of course ruling out the existence of an intrinsic 
double transition doesn't imply that its superconductiv- 
ity is conventional. In particularly, following the idea 
that T C 2 is the intrinsic transition, its strong sample 
dependence (1.76 K for L3-1, 1.71 for Sl-1, 1.68 for Ll-7 
and 1.53 K for L2-2, i.e. a dispersion of 15%) might 
point to unconventional superconductivity. Nevertheless 
this observation must be carefully investigated as the 
dispersion of T c \ is only of 5% for all our samples. 

So we have accumulated evidence for an extrinsic dou- 
ble superconducting transition in PrOs4Sbi2. Thus, in 
addition to account for the two steps visible in magnetic 
or resistivty results, the presence of normal part in the 
sample above T c i would explain the highest temperature 
minima in flux flow resistance reported in reference 2 ^. 

It happens that a definitive conclusion may emerge 
from a quantitative measurement of the specific heat of 
a single transition sample, specially by validating our 
semi-quantitative observation that the single jump at T c 
is higher than the sum of the two jumps in the double 
transition samples. But finding the origin of the extrinsic 
superconducting transition would provide the definitive 
answer. 



B. Origin of the Double superconducting transition 
: discussion and 4 circles X-rays diffraction study 



Our last proof is presented figure [7| showing the ac spe- 
cific heat of sample Sl-1 before and after polishing down 
to 120/^m (then called Sl-lA). As the samples are too 
small to determine their mass, the curves are normalized 
so the entropies in the normal phase match. The data fit 
in all the temperature range from 1.3K to 4K except in 
the double transition. By polishing, a dramatic decrease 
of the ratio of the specific heat jumps occured, going 
from 1.9 to ~1.0. The transition at T C 2 became steeper 
and sharper whereas we aren't even able to distinguish 
any change of slope between T c \ and T C 2 in Sl-lA. So 
polishing the sample Sl-1 clearly tended to remove the 
transition at T c \. 

Finally all these observations are compatible with a 
double transition due to two parts of the sample. The 
effect of polishing presented here as well as reported ini£ 
(we got a single transition sample L1-1A by reducing the 
thickness down to 50/im) suggest that the two parts are 
macroscopically segregated. 

We would like again to draw attention to the problem 
of the quality of the samples and on possible spurious 
analyses of the temperature dependence in the supercon- 
ducting state of the sample with a large single transition. 
Indeed from the specific heat measurements of batch L2, 
it appears that some part of the samples can still be- 
come superconductor at temperature much lower than 
T c i (1.52 K for batch L2) and it would not be surprising 
that it is also the case in the samples with a large single 
superconducting transition. 



As discussed above, the low temperature supercon- 
ducting transition T c i seems to be the intrinsic one. So 
a simple random-impurity-induced pair breaking mech- 
anism cannot be responsible for the appearance of the 
extrinsic transition which occurs at higher temperature 
(T c i) and, whatever the intrinsic transition is at T c \ or 
T C 2, the narrowness of the lowest-temperature jump rules 
out this hypothesis, since such an effect would simultane- 
ously broaden and lower T c . So it is not surprising that 
annealing the sample has no effect on the double transi- 
tioni^.We note that this case is not isolated in the history 
of heavy fermions superconductors. The double transi- 
tion of URu2Si2 was ruled out by Ramirez et alJ^ by 
isolating the lowest T c phase when removing the surface 
of the sample. As for CePtaSi, Kim et al. pointed out 
to a spurious double transition 23 due to a second phase 
of Ce 3 Pt 23 Siii, and T c decreases from 0.75 K to 0.46 K 
with increasing quality of the samplers. Sr2Ru04 — and 
Celrlns 26 exhibit a much higher T c in resistivity mea- 
surements than on the specific heat results. 

Multiple scenario are possible like the existence of an 
impurity phase very similar to PrOs4Sbi2 and supercon- 
ductor at T c i, or the presence of twins which enhance 
T,2L. We test here another hypothetic scenario involving 
the existence of praseodymium vacancies in the structure, 
in relation with the disparity of the quantitative value of 
the specific heat on the Schottky anomaly. Actually in 
the filled skutterudite structure RT4X12, because of a 
weak interaction of the R-atoms with its neighbours (as 
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points out by the large rattling of the R-atoms in the 
X12 cages), some R- vacancies are commonly observed. 
Moreover as discussed in section II, the large dispersion 
of the specific heat value around 2 K could be ascribed to 
varying ratio of Pr-vacancies which may strongly affect 
the Schottky anomaly. Based on these observations, we 
suggest the overall scenario discussed ir>i£, namely the 
transition at T c i appears in the part of the samples with 
only partial occupancy of the Pr site. 

In order to test this hypothesis, we have selected three 
single crystals with different shapes of the superconduct- 
ing transition and we have submitted them to a single 
crystal x-ray diffraction experiment. Fig.[8]shows the su- 
perconducting transition in the specific heat normalized 
to its value in the normal state of these three samples. 
From sample L3-2 to Ll-3, the ratio A(C/T)i/A(C/T) 2 
decreases strongly from 1 to 0.29 reaching zero in sample 
Ll-5. If the hypothesis of Pr vacancies is true, the jump 
at T cl should increase when the vacancies on the Pr sites 
becomes larger. 




FIG. 8: Zoom on the superconducting transition in the spe- 
cific heat normalized to its normal part obtained under mag- 
netic field for the three samples we have selected (for their 
different ratios of the 2 specific heat jumps) for the 4 circles 
X-ray diffraction. Results of 4 circles X-ray diffraction are 
presented table ITT1 

The X-ray investigation was carried out with a No- 
nius KappaCCD diffractometer equipped with graphite 
monochromatized AgKa radiation. After sample align- 
ment, up to 20000 Bragg reflexions were collected to a 
maximum sin(9/A of 1.15 leading to a very high redun- 
dancy. After extraction of the intensities using the Eval- 
CCD software^, a numerical absorption correction was 
applied using the crystal shape. The structure refine- 
ment was carried out using the Jana2000 software^. An 
isotropic extinction correction (type I, Lorentzian distri- 
bution) was applied and all atoms were given anisotropic 



S ample 


TOO 

LiO-Z 


T1 o 
Lil-o 


T 1 ^ 

J_j1-0 


Cell parameter (A) 








ripp ( Pr 1 

UL.L . 1111 


0.93(1) 


0.89(1) 






u.uo i vyo j 


u .uoocJy i j 




xj iso (Os) (A 2 ) 


0.0046(1) 


0.00402(5) 


0.00480(3) 


u lso (Sb) (A 2 ) 


0.0064(1) 


0.00578(6) 


0.00660(5) 


x(Sb)(i) 


0.15608(5) 


0.15613(3) 


0.15608(2) 


z(Sb)(i) 


0.34040(5) 


0.34036(3) 


0.34031(2) 


Gof 


2.24 


1.95 


2.33 



TABLE II: Structural parameters and refinement agreement 
factors for three crystals. The position of Pr atoms is setting 
at (0,0,0). Uiso : isotropic thermal parameters. The occu- 
pancy (occ.) of Os and Sb was set to 1. Gof is the goodness 
of fit. 



atomic displacement parameters (a.d.p.). Finally, since 
an anomalously large a.d.p. was observed for the Pr atom 
(about 0.04 A 2 ) and to test the vacancy ratio on the Pr 
sites, its occupancy factor was also let to vary. This 
systematically led to a slight decrease of the Pr atom oc- 
cupancy (from 1 to 0.97 for sample Ll-5) and the a.d.p. 
remained practically the same. The agreement factors 
were improved, though only slightly : for crystal Ll-3 
having the lowest refined Pr occupancy (0.89), the good- 
ness of the fit decreased from 2.05 to 1.95 by letting the Pr 
occupancy parameter free. Table [TT] reports the param- 
eters obtained for the three single crystals when the Pr 
occupancy is refined and the Pr position is set at (0,0,0). 

These result confirm that the Pr atom is in a strongly 
disordered position as reported in&2ii. Moreover, they 
seem to indicate the possibility of a nearly but not com- 
pletely filled Pr site in the filled skutterudite structure for 
PrOs4Sbi2- In the interests of caution, we should point 
out that the large rattling of the Pr in the Sb-cages at 
300 K may reduce the accuracy of this measurement, and 
the values may depend slightly on the refinement model. 
However the occupancy factors reported in table [TTJ at 
least support the trend of varying Pr occupancy in these 
three samples and the Pr occupancy is the highest in 
sample Ll-5 so with a single sharp transition. However a 
smaller level of vacancies in sample L3-2 than in sample 
Ll-3 is not relevant with the above-mentioned hypothesis 
since sample L3-2 has a larger ratio A(C /T)i/ A(C /T)%. 
We think that further diffraction measurements at low 
temperature are required to thoroughly characterize the 
dynamic and static disorder at the Pr site, and possi- 
bly prove or disprove the hypothesis Pr vacancies as the 
origin of the double transition. 



C. Superconducting phase diagram 

We have followed the superconducting transition of the 
single transition sample Ll-5 under magnetic field by re- 
sistivity measurements. In figure 9(b) we report T c ver- 
sus H, determined by onset criterion. The transition 
remains very sharp (12 mK at 1.2 T and less than 30 mT 
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at 400 mK) pointing out again to the high quality of this 
sample. The T C (H) line matches with T c2 {H) published 
in^ as shown in figure 9(a) The small positive curva- 



ture at low magnetic field is even more clearly visible (Cf. 
insert of figure 9(a)). 
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FIG. 9: (a) Superconducting phase diagram under magnetic 
field of the single transition sample Ll-5 and the sample with 
a double transition from^i. The insert presents a zoom on 
the small positive curvature of both the single transition sam- 
ple Ll-5 and T C 2(H) of the sample with a double transition 
frorn^. Clearly the superconducting transition of Ll-5 follows 
the lowest-temperature one of the double transition sample. 
The small positive curvature at low field is even better re- 
solved, (b) Superconducting phase diagram under magnetic 
field of a sample with a single sharp transition Ll-5 obtained 
by resistivity measurements. The error barres indicate the 
witdh of the transition. The data were fitted by a two-band 
model like in refi^. The new set of parameters are described 
in the text. The dashed dotted line shows the fit without 
Pauli limitation (with g=0)^i. 



This makes it even clearer that multigap effects are 
disconnected from the question of the double transition. 
It also backs up the conclusions of the recent thermal 



conductivity measurements on a high quality-single su- 
perconducting transition sample, Ll-lA: multigap effects 
have been confirmed^, with a very low field scale asso- 
ciated with the light carrier/small gap band, of the same 
value as found in a previous inhomogeneous (wide spe- 
cific heat transition) sample^. Moreover, the fact that 
the small positive curvature close to T c is found also in 
homogeneous samples and with the similar amplitude is 
a definite proof, beyond the reproducibility of the mea- 
surements^, that it is not connected to sample inhomo- 
geneities. 

It also brings information on the inter band coupling 
strength different from those of the thermal conductiv- 
ity experiments. For example, the fit proposed ini^ for 
H C 2(T) would also apply to these new measurements, as 
the data simply scale with T c . But the set of inter and in- 
tra band couplings (Ay ) proposed in this first work, was 
based on the simplest hypothesis that \j is proportional 
to the density of states of band j, so that An = A21 and 
A12 = A22- In such a case, a simple calculation of the 
two gaps in a weak-coupling scheme shows that they are 
equal. In order to be consistent with the thermal con- 
ductivity results, which find a factor three between the 
small and large gapi^, one needs to introduce a differ- 
ence between An and A21. The size and position of the 
curvature on H& then still impose a very small value of 
A12 (we take still for simplicity (A12 = A22). 

So, instead of the set of parameters : An = A21 = 1, 
A12 = A22 = 0.04 proposed iniS, we propose the new set 
: An = 1, A21 = 0.2, A12 = A 2 2 = 0.07 and g— 2, w hich 
yields a fit of the same high quality (see fig. |9(b)[ ), but 
yields also the good values for the gap as deduced from 
thermal conductivity measurements. Again, it is only 
the ratio of the Ay which matters, the value An = 1 
being arbitrarily fixediS. The factor 5 between An and 
A21 is essentially due to the coupling strength, meaning 
that intra band coupling in the band with heavy effec- 
tive masses (having f character) is much stronger than 
inter band coupling from this band to the band with a 
small mass (weak f character) . Of course, A12 and A22 are 
strongly reduced by density of states effects, but the gen- 
eral trend which emerges from the new set of A,j imposed 
by the combination of thermal conductivity 19 and H C 2 re- 
sults is that multiband effects in PrOs4Sbi2 are coming 
from the difference in the f character of the bands both 
through density of states and pairing mechanism effects. 
This conclusion is quite robust as it relies on measure- 
ments independent of the sample homogeneity and the 
number of transitions. 



IV. CONCLUSION 

Although general characterizations point out to a re- 
current double superconducting transition in PrOs4Sbi2 
appearing in the samples with the best RRR and high 
specific heat jump at the superconducting transition, so 
in good samples, a study of many samples specially small 
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ones (with a typical size of 100 /im) shows that its oc- 
currence could be a phenomenon related to an inhonao- 
geneous effect rather than to fundamental microscopic 
mechanisms. More precisely we think the lowest temper- 
ature transition T c2 is the intrinsic one. Based on our 
4 circles X-ray diffraction results, we conclude that Pr 
vacancies are certainly present in the samples and with 
various percentages which might explain the discrepancy 
between quantitative specific heat results. But estab- 
lishing a clear relationship between Pr vacancies and the 
occurence and magnitude of the jump in specific heat at 
T c i needs further studies, specially at low temperature. 

Finally, the superconducting phase diagram of a single 
transition sample was determined and fitted with a two- 
band model, ft appears in connection with thermal con- 
ductivity results that the multiband effects in PrOs4Sbi 2 
come from the difference in the f character of the bands 
both through density of states and pairing mechanism 
effects. 
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